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LONG-TERM  GOALS 

The  goal  of  this  project  is  to  examine  the  role  of  sea-ice  and  atmospheric  interactions  in  the  retreat  of 
the  SIZ.  As  sea  ice  retreats  further,  changes  in  lower  atmospheric  temperature,  humidity,  winds,  and 
clouds  are  likely  to  result  from  changed  sea  ice  concentrations  and  ocean  temperatures.  These  changes 
in  turn  will  affect  the  evolution  of  the  SIZ.  An  appropriate  representation  of  this  feedback  loop  in 
models  is  critical  if  we  want  to  advance  prediction  skill  in  the  SIZ.  To  do  so,  we  will  conduct  a 
combination  of  targeted  measurements  and  modeling  experiments  as  part  of  the  atmospheric 
component  of  the  Seasonal  Ice  Zone  Reconnaissance  Survey  project  (SIZRS). 

OBJECTIVES 


•  Determine  the  role  of  changing  atmospheric  properties  in  modifying  the  evolution  of  the  SIZ  in 
the  Beaufort  and  Chukchi  Seas  from  spring  through  fall. 

•  Determine  how  changes  in  sea  ice  and  sea  surface  conditions  in  the  SIZ  affect  changes  in  cloud 
properties  and  cover. 

•  Determine  the  role  additional  atmospheric  profile  observations  may  play  in  improving  the 
quality  of  weather  forecasts  and  ice  predictions  for  the  SIZ  of  the  Beaufort  and  Chukchi  Seas. 

•  Adapt  a  low  cost,  expendable,  air-deployed  micro-aircraft  to  obtain  temperature  and  humidity 
profiles  and  cloud  top  and  base  heights 
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APPROACH 


To  achieve  these  long-term  objectives  we  will  conduct  observation  and  model  experiments.  The  SISZR 
project  is  an  integrated  observation  and  modeling  program  aimed  at  understanding  the  interplay  of 
atmosphere,  ice,  and  ocean  in  the  SIZ  of  the  Beaufort  and  Chukchi  seas  (BCSIZ).  It  will  take 
advantage  of  routine  Coast  Guard  C-130  domain  awareness  missions  that  take  place  at  two-weekly 
intervals  from  March  through  November.  As  part  of  the  atmospheric  observation  component  of 
SIZRS,  this  project  will  deploy  dropsondes  during  SIZRS  flights  planned  at  least  monthly  from  April 
through  October  to  obtain  atmospheric  profdes  of  temperature,  humidity,  and  winds.  Cloud  top  heights 
will  be  retrieved  using  infrared  imagers  and  LIDAR  deployed  aboard  the  SIZRS  aircraft  by  other 
SIZRS  projects.  Sea  surface  temperatures,  ice  concentrations,  and  floe  size  distributions  will  be 
measured  by  other  components  of  the  SIZRS  project  as  well.  Our  atmospheric  observations  will  be 
examined  in  the  context  of  varying  surface  conditions  (sea  ice  concentration,  ice  thickness,  and  SST) 
to  increase  our  understanding  of  atmosphere-ice-ocean  interactions  and  to  initialize,  validate,  and 
improve  our  meso-scale  atmospheric  model.  Seasonally  changing  surface  conditions  are  expected  to 
provide  a  present  day  analog  for  expected  future  ice  retreat.  Ultimately  we  hope  to  transmit  dropsonde 
observations  to  forecast  centers  for  potential  assimilation  into  analyses  and  forecast  products.  In 
addition,  we  will  contribute  to  technology  development  by  adapting  and  deploying  a  new  generation  of 
truly  expendable  (<$700)  micro-aerial  vehicles  (Data  Hawk  and  SmartSonde)  designed  to  obtain 
detailed  high-vertical-resolution  temperature,  humidity  and  wind  profdes  and  cloud  layering 
information  that  cannot  be  obtained  with  traditional  dropsondes.  Our  vision  is  that  these  vehicles  will 
deliver  new,  inexpensive  measurement  capabilities  for  research  and  operational  purposes  in  the  data 
sparse  region  of  the  BCSIZ  as  well  as  other  regions  of  the  globe.  Ron  Lindsay  (UW)  and  Zheng  Liu,  a 
recenty  hired  post-doc  contributes  WRF  modeling  expertise,  Jinlun  Zhang  (UW)  conducts  sea  ice 
model  experiments.  Dale  Lawrence  and  Jim  Maslank  (CU)  contribute  to  AUV  expertise. 

WORK  COMPLETED 

Observations: 

•  We  adapted  a  commercial  GPS-based  radiosonde  system  to  operate  in  a  dropsonde  mode  that  can 
be  launched  from  aircraft. 

•  We  completed  the  Aircraft  Configuration  Control  Board  (ACCB)  process  including  safety  of  flight 
test  (SOFT)  and  obtained  approval  for  deploying  dropsondes  during  ADA  flights. 

•  We  conducted  successful  deployments  of  dropsondes  during  ADA  missions  in  June,  July  and 
August  with  a  total  of  26  profiles  over  the  course  of  6  days. 

•  We  identified  and  solve  several  problems  associated  with  the  dropsonde  deployment  method.  We 
experimented  with  different  devices  used  to  slow  the  descent  of  the  dropsondes  and  arrived  at  a 
solution  that  provided  desireable  descent  rates  and  allowed  deployment  within  the  constraints  of 
ADA  missions. 

•  We  intercompared  satellite  retrievals  of  cloud  fraction  and  cloud  top  height  with  observer  estimates 
from  C-130  cockpit. 

•  We  assessed  the  potential  of  simulatenous  CloudSat/Calipso  retrievals  with  dropsondes  taking  into 
satellite  sampling  patterns  and  spatial  and  temporal  variability.  This  information  will  help  guide 
future  mission  planning. 
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We  worked  towards  IR-camera  based  solutions  to  obtaining  cloud  top  height  as  an  alternate  to  the 
CUPLIS-X  system  (PI,  Mark  Tchudi,  University  of  Colorado)  for  which  approval  is  still  pending. 
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Fig.  1.  Terra  MODIS  satellite-derived  cloud  fraction  and  cloud  top  pressure  over  the  WRF 
model  domain  (30  km  resolution)  around  2200Z  during  the  August  16,  2013  SIZRS  mission.  The 
boundaries  of  the  10  km  resolution  nest  domain  is  shown  in  thick  black  lines.  The  sea  ice 
concentration  contours  are  shown  in  thin  black  lines.  The  location  of  dropsondes  deployment 

are  shown  in  green  dots. 


Modeling: 

•  A  post-doctoral  research  associate,  Zheng  Liu,  was  hired  to  conduct  WRF  model  experiments 

•  We  conducted  Weather  Research  and  Forecast  model  simulations  (WRF)  the  summer  of  2013  and 
conducted  experiments  using  different  resolutions,  initalization  schemes  and  parameterizations. 

•  We  intercompared  different  model  runs  with  SIZRS  dropsonde  observations  and  the  NCEP  Global 
Forcecast  System  (GFS)  and  ERA-Interim  reanalysis  data.  We  examined  differences  in 
temperature  and  humidity  profiles,  horizontal  wind  speed  and  the  role  of  clouds  in  the  differences 
between  models  runs  and  observations. 
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Advanced  Observation  Playforms  (DataHawk,  SmartSonde) 

•  Work  on  a  SmartSonde  design  to  obtain  detailed  atmospheric  parameters  and  cloud  top  and  base 
and  can  be  launched  from  C- 130  is  progressing.  A  more  detailed  report  is  submitted  separately  by 
Co-Investigators  Lawrence  and  Maslanik. 

RESULTS 

Near  surface  temperature  inversion 

The  June  SIZRS  temperature  profiles  are  characterized  by  a  strong  near  surface  temperature  inversion 
as  is  shown  in  Fig.  2.  The  temperature  profiles  from  the  WRF  simulations  and  the  NCEP  Global 
Forecast  System  (GFS)  and  ECMWF  ERA  Interim  (ERA-I)  reanalysis  reproduced  this  feature.  We 
followed  the  Kahl  (1990)  method  to  identify  the  inversion  layer  and  use  the  temperature  change  across 
the  inversion  (inversion  strength)  to  assess  the  modeled  inversion. 
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Fig.  2.  Comparison  of  mean  and  R.M.S.E  of  temperature,  relative  humidity  and  horizontal  wind 
speed  for  the  June  2013  SIZRS  campaign.  The  WRF  simulated  profiles  are  shown  in  blue,  with 
solid  lines  for  the  parent  domain  and  dashed  lines  for  the  nested  high  resolution  model  domain.  The 
dropsonde  measurements  (snd),  the  GFS  and  the  ERA-Interim  (ERAI)  reanalysis  data  are  shown  in 

red,  black  and  green  respectively. 
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The  inversion  strength,  on  average,  is  overestimated  in  WRF  simulations  by  1 . 1  K  but  underestimated 
by  the  GFS  and  ERA-I  reanalysis  by  2.1  K  and  1.4  K  respectively.  The  root  mean  square  (RMSE) 
error  of  WRF  simulated  inversion  strength  is  1.7  K,  also  smaller  than  the  reanalysis  (3K),  which 
indicates  additional  skill  by  the  local  WRF  implementation.  The  overestimated  inversion  strength  by 
WRF  is  related  to  a  very  shallow  layer  of  cloud  close  to  surface,  lower  than  100  m.  It  raises  the  base 
height  of  inversion  to  its  cloud  top  level  from  the  surface  and  introduces  a  cold  bias  to  the  inversion 
base  temperature  by  radiative  cooling  up  to  3.5  K. 

Bias  in  reanalysis  datasets 

The  August  SIZRS  mean  temperature  and  relative  humidity  profiles  suggest  a  low  level  cloud  layer 
embedded  in  a  deeper  weak  inversion  layer  (Fig.  3).  However  this  feature  is  not  captured  by  the  GFS 
reanalysis,  which  has  a  systematic  moist  bias  over  25%  RH  above  cloud  to  780  hPa.  Note  that  WRF 
simulated  RH  profiles  has  a  slightly  smaller  bias  and  RMSE  than  the  GFS  reanalysis,  which  suggests 
that  although  WRF  has  the  potential  to  improve  the  forecast,  the  bias  in  the  forcing  dataset  is  still 
critical. 


Fig.  3.  Comparison  of  mean  and  R.M.S.E  of  temperature,  relative  humidity  and  horizontal  wind 
speed  for  the  August  2013  SIZRS  campaign.  The  WRF  simulated  profiles  are  shown  in  blue,  with 
solid  lines  for  the  parent  domain  and  dashed  lines  for  the  high  resolution  nest.  The  dropsonde 
measurements  (snd)  and  the  GFS  reanalysis  data  are  shown  in  red  and  black  respectively. 
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Another  systematic  bias  in  reanalysis  is  the  horizontal  winds.  The  wind  speeds  are  biased  low  for  most 
levels,  especially  the  low  level  jet  (Fig  2.  and  Fig.  3).  The  bias  is  also  evident  in  the  meridional  wind 
speed,  which  is  persistent  from  surface  up  to  700  hPa  (Fig.  4).  Since  horizontal  winds  are  nudged  at  all 
levels  in  WRF  simulations,  it  is  not  likely  for  WRF  to  drift  far  away  from  the  forcing. 


Fig.  4.  Comparison  of  mean  |t/|  and  \V\  for  June  (a)  and  August  (b)  SIZRS  campaign.  The  WRF 
simulated  profiles  are  shown  in  blue,  with  solid  lines  for  the  parent  domain  and  dashed  lines  for  the 
high  resolution  nest.  The  dropsonde  measurements  (snd),  the  GFS  and  the  ERA-Interim  (ERAI) 
reanalysis  data  are  shown  in  red,  black  and  green  respectively. 

Model  resolution  and  parameterizations 

Our  WRF  simulations  have  a  parent  domain  with  a  horizontal  resolution  of  30  km  and  a  nest  domain 
with  a  horizontal  resolution  of  10  km  (Fig.  1).  The  near  surface  temperature  and  moisture  biases  are 
slightly  reduced  for  the  June  SIZRS  sounding,  using  the  higher  horizontal  resolution.  For  the  August 
SIZRS  soundings,  the  higher  horizontal  resolution  shows  little  advantage.  A  set  of  WRF  simulations 
with  doubled  number  of  vertical  levels  are  examined.  The  additional  levels  have  very  limited 
improvements  compared  to  the  low  resolution  runs  (not  shown). 

We  experimented  with  different  choices  of  parameterizations  in  WRF,  e.g.,  Goddard  microphysics  vs. 
Morrison  microphysics,  MYJ  vs.  YSU  boundary  layer  schemes  and  Grell-Devenyi  vs.  new  Simplified- 
Arakawa-Schubert  cumulus  parameterization.  We  found  that  the  inter-model  variability  is  much 
smaller  than  the  model  bias  with  respect  to  the  dropsonde  observations.  For  example,  the  June  SIZRS 
temperature  profiles,  the  maximum  inter-model  standard  deviation  is  1 .2  K,  compared  to  the  maximum 
model  bias  standard  deviation  of  5.3  K.  Further  detailed  analysis  is  required  to  identify  an  optimal  set 
of  parameters. 

IMPACT/APPLICATIONS 

The  results  point  the  way  for  future  improvement  in  model  simulations.  Biases  in  the  inversion 
strength  simulated  by  models  and  analysis  have  potentially  important  implication  for  the  transfer  of 
energy  from  the  atmosphere  to  the  sea  ice  during  the  melt  season. 
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The  characterization  of  biases  in  model  and  analysis  wind  speeds  is  important  for  sea  ice  model 
development.  Sea  ice  models,  driven  by  winds  from  atmospheric  analysis  are  senstitive  to  errors  in 
surface  winds. 
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